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ABSTRACT 


Simultaneous measurements of water surface fluctuations 
mieenhorizeleal water particle velocities in a line perpen- 
dicular to the direction of wave approach extending across 
mee suri Zone were taken in varying surf conditions at 
mo locations. tlhe Spectral velocities calculated using 
linear theory as a transfer function underestimated measured 
values by 79-86% at the peak of the spectrum. The coherence 
values were generally low indicating non-linear and turbulent 
conditions. Strong harmonics in the spectra of the waves 
pia Water Darticle velocities further suggest a non-linear 
System. The eenee eal phases computed using linear 
theory did not accurately predict the observed phases. In 
general breaking waves can be characterized as a strongly 
non-linear wave phenomenon. Measured frequency distributions 
were compared with both Gaussian and Gram-Charlier distribu- 
meoins Dy USing the chi-square goodness-of-fit test. Quali- 
macively, the Gram-Charlier distribution gave the better fit 


to the flow velocity data. 
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Pee PRODUCT TON 


Wave theories developed for deep water waves can be 
applied with some degree of certainty, and can be tested 
in laboratory and field situations. These theories general- 
meerold until such time that the shoaling process begins, 
and with somewhat lesser accuracy, throughout the shoaling 
eocess up to the point of near breaking. At the breaker 
Pemamt, however, there 15 a transition from ordered to 
Perpulent motion and the description of wave kinematics 
becomes more difficult. The most forthright approach to 
mae problem of describing the kinematics of surf zone 
breaking waves is through direct measurements. Advances in 
micerument design have led to simple, sturdy, measuring 
devices with rapid response time which can senSe the small 
scale as well as the large scale features. 

The study of the kinematics of breaking waves in the 
Surt zone has progressed slowly due both to the problems 
pEmecuntercd in making direct field measurements and the 
mmerroulty in modeling the surf zone in the laboratory. 
Breaking waves were measured in the laboratory in early 
work by Iverson (1953) using photographic techniques. 

Me length of the channel limited the wave type to plunging 
and surging breakers. Adeyemo (1970) made similar labora- 
tory measurements using hydrogen bubble and photographic 


methods. Gaughan and Komar (1975) applied the theory of 
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Wave propagation in water of gradually varying depth, as 
developed by Biesel, to determine the dependence of breaker 
type on the beach slope tangent and the deep water wave 
Seccpness. 

Inman and Nasu (1956) made field measurements of water 
particle velocity by measuring the drag force under the 
Wave in order to infer particle velocities. Miller and 
Ziegler (1964) used both acoustic and electromagnetic 
current meters to determine the particle motion in the 
surf zone. Walker (1969) made studies using propeller-type 
flow meters. Wood (1973) measured waves and currents in 
the surf zone using movies of dye movement and capacitance 
Wave gauges. Fuhrboter and Busching (1974) utilized a two- 
component current meter and two pressure-type wave meters 
to measure simultaneous orbital velocities and water levels. 
Mim@ernton (1968), Steer (1972), Thornton and Richardson 
feo) and Bub (1974) used pressure meters, capacitance 
Wave gauges and electromagnetic current meters to measure 
emacs protiles and particle velocities; the work 
meesented here is an extension of these studies. 

IMewooVeCtEIVe Of this rescarch 1S to study the kine- 
Weietes Of Water particles in breaking waves within the 
Suri zone. Simultaneous measurements were made of the 
instantaneous wave profile and the horizontal water 
Mametec le velocities at fixed locations in the surf zone. 


feimates of the probability density functions and power 


tal. 








Spectra of the wave heights and particle velocities are 
made. The applicability of using linear wave theory as 
peeopectral transfer function in computing velocity spectral 
components from the power spectrum of the waves is 
measured. The computed velocity spectra are compared with 
actually measured velocity spectra. Theoretical phase 
Spectra between wave gauges are calculated and compared 


to measured phase spectra. 


a2 








II. MEASUREMENTS 


PweeeEXPERIMENT SITES 

the experimental, sites are in the vicinity of Monterey, 
California. The beaches here are some of the first inten- 
Seavely studied to gain an understanding of amphibious war- 
fare techniques and are described by Bascom (1964). It was 
meoired tO Measure the various types of breaking waves in 
Sie Study including spilling, plunging and surging breakers. 
The manner in which waves break depends very much on the 
characteristics of the beach and near-shore bottom slope 
fj oole 1). Plunging and spilling breakers occur most 
frequently at the Del Monte Beach site within Monterey Bay. 
The waves at this location are generally severely direc- 
muemaliy tiltered and refracted by the geometry of the bay 
and impinge almost perpendicular to the shore with a 
resulting simplification in the wave description. A second 
experimental site was Carmel River Beach, five miles to 
the south, where the beach is very steep and very often 
has surging type breakers. Again here, the beach is within 
an embayment and the waves impinge almost perpendicular to 
mreeshore. A typical beach profile and instrument loca- 
tion for Del Monte Beach is shown in Figure 1. Del Monte 
Beach profiles for 4 - 10 March 1975 are given in 
Appendix A. Figure 2 shows the beach profile and instrument 


Peceiion at Carmel River Beach on 29 May 1975. A beach 


eS 











SLoqzew 7Z 02 dn 


Spliees cawair 


SUT[TTTds pue SuUTSins 


mn Otp 


tedsp[oy z31eNdH 


C2 etine tl iene: 2 


SL6L ACW 62 


ydeog LOATY Towase) 


19}0w [T 03 dq 


Spuodas OT OF GOT 


SUTTTTds pue sutsunNtTg 


n 00S 


sedsp[To}F Z}321eNDH 


VY xtpuoddy pue [ oansty 


SZ6T YItLEN OT - P 


yoeog o2UOW T9qG 


SJILSIUSLIVAVAD aAVM GNV HOVSd 


TYSTOH OALH 


potiog oACA 


odd] 1LoOyeoig 


LoZowerqg ueIpsy pues 


addy, pues 


aedo[s yoreg 


93eq 


UOTIEDOT 


eae anya: 


14 








SZ6T YoOteW OT - PV 
yoeoeg o2UOW [T9q 28 UOTeEIOT YUOUNASUT pUe sTTFoLIg yoeeg [edtdA], “TL AUNnoIA 
SUILIN NI SJONVLISIA 
. 09- 0S- Op- 0c - 0Z- OT- 0 
fe 
T° OL NOILVUFIOVXA IWOILYUAA 
(IA HIGVL) AdOIS HOVAG 


TH (MTIN) 
UILAN 


MOTS = | 








SUN 


Lay) 





¢# HONVO HAVM 


C# FONVI FAVM 





0 


OOT+ 


UOGs 


O0¢+ 


Sameer Ni ONO IVAW 1S 


1 





SL6L AGN 6Z 
Yoeog LIATY LTouwle) a is Ss Ort. Og Pou es Sy pue SL CFOLd yoeog aes TUnOI[A 


SUSLHN NI SONVLSIC 
No O.T- 0 OT+ O¢+ OC+ 





TO NOL yaa ooN a oy oun ay 
Eo CeeetdOrs Kova 


li 
YALAN 
MO'TA 





2 x + 





Ay (Z@# «TOYO SAM 
IO WOLLOG) 


[T# 4ONVD SAVM— 


_/—2# s9onvo sAvM 


OOT- 


OOT+ 


OU 


Oey 


SUALAWL INA NI 


16 





Meorrle taken the following day was identical to 


Eegure 2. 


B,. INSTRUMENTATION 

iainy problems were €mcountered in the initial attempt 
imjeeebruary 1975 to measure the desired parameters. The 
first wave gauges were manufactured from 1.0 cm diameter 
Serapped RG-11 coaxial cable as was done by Bub (1974). 
Large masses of kelp torn loose by high waves, flotsom 
and the pounding surf easily broke these cables. Scouring 
around the base-plates of the instrument towers and the 
blade-type anchors caused the towers to topple. The con- 
struction of the wave gauges was modified to a three- 
eaehths inch outside diameter stainless steel rod, covered 
With a one-sixteenth inch wall thickness teflon tubing 
fitted tightly over the rod; these wave gauges withstood 
the severe forces in the surf well. Steel pipe extensions 
each 0.6 m in length were added under the tower base plates 
and resulted in stable platforms. An extension was built 
for the top of the anchors; with this device, the anchors 
mere Siledge-hammered deeper into the beach so that the top 
of each anchor was approximately 15 cm below the sand. The 
extension was then removed and scouring was minimized. 

Measurements at Del Monte Beach were made using two 
mmeetromagnetic flow meters and three capacitance wave 
gauges. One flow meter and two wave gauges were used at 


Carmel River Beach. The instruments were mounted on the 


ily 








towers Within the surf zone. All equipment was calibrated 
imeecthe laboratory prior to the experiments, 

MitenmclCctrenics package for each of the 2.5 m 
Capacitance Wave gauges was constructed using the design 
of McGoldrick (1969). The capacitance wave gauge operates 
@amene principle that a change of capacitor plate dimensions 
@nemees the circuitry voltage output. The rod and the 
sea water act as plates and the teflon insulation as the 
@memectric. As the wave elevation fluctuates, the capaci- 
memee Of the circuit changes and the voltage output responds 
[Memearly. These water surface fluctuations were sensed by 
fmemansisStorized circuit operated by self-contained bat- 
teries. The electronics package was housed in a water- 
tight brass case which was mounted on the tower. This 
eeaaeled the connecting leads to be less than 30 cm, thereby 
Minimizing Wire-to-wire capacitance. Accuracy was estimated 
pombe +0.005 m. The calibration plots are shown in 
Appendix B. Minor modifications to the wave gauges 
immediately following the Del Monte Beach measurements 
@eewtnie cause for the different calibration graphs in 
mopendix B. 

The two flow meters were Marsh-McBirney Model 721 
meeetromagnetic Current Meters. The flow meter operation 
is based on Faraday's principle of electromagnetic induction. 
Each probe measures water velocity in two orthogonal 


frections. Fhe flow meters were calibrated with an 


18 





Seertliiating platfiorm attached to an eccentric arm driven 
Dy a variable speed motor using the method of Thornton 
and Krapohl (1974). Measurement accuracy was determined 
to be +0.005 m/sec during calibration. The flow meter 
Subapration plots are shown in Appendix C. 

Mnethe field the instruments were attached to 6.3 cm 
eeraiae diameter steel pipes which were 3.6 m high with a 
1.0 m base plate and the 0.6 m steel pipe extensions. The 
ewers were placed on a line perpendicular to the predomi- 
nant wave direction and were erected during low tide when 
the beach was easily accessible. The measurements were 
miemeconducted at high tide. The tidal range in the 
Monterey area 1s typically two meters. The towers were 
elperied Dy Steel puy wires fastened to the anchors. 
mBeveral types of anchors were tried and the blade anchor 
worked quite satisfactorily in the relatively coarse sand. 
The wave gauges were attached to the towers with locally 
fabricated micarta clamps. A1.6cm Steel stiffening rod 
was added to the arrangement in order to reduce wave gauge 
vibration following breaker impact. The flow meters were 
positioned directly under the wave gauges with the axes 
@ieponed horizontally (+X, shoreward) and vertically (+Y, 
upward). A carpenter's level was used to establish axis 
alignment with an estimated error of +2 degrees. A typical 
tower and sensor arrangement is shown in Figure 3. 

All signals were cabled ashore and recorded on a Vidar 


Sevoratron, 5Z2-Channel! digital data acquisition system. 
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A Sangamo model 3500 14-channel FM tape recorder was 
utilized as a secondary recording system. A Brush 8-channel 
strip chart recorder was used to monitor the instrumentation 
Periormance during recording and as a means to select the 


appropriate data sections to be analyzed. 
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IIIf. ANALYSIS OF DATA 


Record lengths of 30 minutes from each data set were 
analyzed. The details of the analysis system including 
Seer oration factors are given in Appendices D and E. A 
mean value was calculated for all data sets and the data 
was linearly detrended to exclude the rise and fall of the 
tide. The variance, standard deviation and average period 
were calculated. The average period was determined by 
calculating the time between zero upcrossings. However, 
perturbations caused by secondary gravity waves "riding" 
the primary gravity waves, capillary waves and minor instru- 
ment noises increased zero upcrossSing occurrences and 
caused the calculated average periods to be lower than 
ferematiy observed. Probability density functions for each 
data set were calculated and graphically compared with 
Samgestan and Gram-Charlier distributions. 

For each signal an auto-covariance function was calcu- 
lated and smoothed with a Parzen window. A Fourier trans- 
form was then applied to the smoothed auto-covariance 
function and the power spectrum determined. A cross- 
covariance function between data setsS was computed, smoothed 
With a Parzen window, a Fourier transform applied, and the 
Suess-Spectrum computed. The coherence and phase were then 
calculated. To compare waves and flow velocity, the wave 


Meorile Spectrum was converted to a theoretical velocity 


ae 








Seectrum £LOr comparison with the measured flow velocity 
meectrum. the two power spectra, the coherence and the 
PpNase were then plotted. In comparing two wave profiles, 
m@eoeretical celerities were calculated using linear theory 
and its shallow water approximation which were then con- 
verted to phases and plotted for comparison with the 
Weaaured phase difference. 

Aliasing testing was performed on all spectra. 
Nyquist frequencies of 1.25 and 0.98 Hz were determined 
optimal for Del Monte Beach and Carmel River Beach data 
Sets, respectively. fhe frequency bands of 0.0 to 1.25 
and 0.98 Hz include the region of gravity waves (0.033 
omen Hz). 

The maximum lag time in calculating the covariance 
fmetyons WaS taken as five percent of the record giving 
meeoeectral bandwidth resolution of 0.00556 Hz and results 
in 40 degrees of freedom for each spectral estimate. The 
80 percent confidence limits for 40 degrees of freedom 
USing a chi-square distribution are found to be between 


teeeorand 1.30 of the measured power spectral estimates. 
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TO SIRI SUI ie 


pee QUALITATIVE DESCRIPTION 

A number of similarities of wave form can be observed 
fmommevariOus types Of breakers occurring on different 
Peeaenes around the world. Figures 4 and 5 are typical 
analog records of waves and horizontal velocities beneath 
the waves obtained from Del Monte and Carmel River Beaches, 
Mmeepectively. In general, there 1S a quick drawdown of 
Pemeer just betore the breaking wave arrives, followed by 
feweccD, almost vertical leading edge, and a sloping pro- 
ies cOward the trailing edge. On the trailing edge of 
the wave secondary waves are often noted. These are 
harmonics of the primary wave frequency and are indicative 
of very non-linear waves. At Del Monte Beach the waves 
break as plunging and spilling breakers giving a generally 
Saw-toothed shape; the plunging or spilling occurs rapidly 
at the crest and moves down the wave. The surging breakers 
occurring at Carmel River Beach rise and fall off more 


gradually. 


meee PROBABILITY DENSITY FUNCTIONS 
fee Gauscstan and) Gram-Chariver Frequency Distributions 
Frequency distributions were computed for each 
measured quantity for all runs and compared with Gaussian 


and Gram-Charlier distributions. The chi-square 
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goodness-of-fit test was computed using both theoretical 
distributions. The variance, standard deviation, skewness 
and kurtosis were also computed for all data sets and are 
Summarized in Table II. The closer the chi-square fit 
Meweaiieter 1S to zero, the better the fit. 

In deep water,the sea surface has approximately a 
GauSsian distribution which results in a skewness of zero 
paeeeeeKUTtosSis Of 3. However, non-linearities introduce 
Skewness and kurtosis values that deviate from Gaussian 
and result in a distribution more closely approximated by 
poe Gram-Charlier. These non-linearities are evident in 
the frequency distribution of Wave Gauge #2 on 4 March as 
Smewneinm Figure 6, and result in bi-modal or tri-modal 
Geer rhutions. Figure 7 is the graph of distributions 
computed for the particle velocity at Flow Meter #2 on the 
Same date. The latter has a uni-modal distribution. 
Pyeemaryx F contains distributions computed for all other 
data sets. 

nae Gastributions for all flow velocities are either 
uni-modal or have less pronounced secondary peaks compared 
Othe distributions of wave heights. As noted in Table II, 
this results in nearly equal fit parameters for both Gaussian 
oma@eeGcram-Charlier distributions, although the latter values 
meuemciialler, and therefore sive the better fit in ten of 
eleven cases. The exception is Flow Meter #2 on 8 March. 
Pagiutionally, the fit parameters for velocity components 


are much closer to zero than are the values computed for 
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wave components. This is to be expected since the flow 
Merers do not experience the surface irregularities to 

the same degree as do the wave gauges. In seven of 
merrteen fit parameter calculations for wave components, 
the Gaussian distribution provided the better fit. This 
1S a result of using a truncated form of the Gram-Charlier 


mecalculating the probability density function, P 


Gee 
given by 
mi m4 ms 

Pac (8) - P.(é) [1.0 ee H. + aq Hy “5 a5 He era } (1) 
where 

Po = Gaussian probability density function, 

mz = Skewness of the data record, 

mM, = Gicmmunrtosas Minus 3, 

em = Hermite polynomials of degree n. 


Higher order moments could be included in equation (1) 
to possibly improve the fit, but would necessitate analyzing 
longer record lengths in order to maintain confidence. 
Eeemronary problems would be encountered in the longer 
fmeeoras, this producing other errors. Hence, the poor 
meets a2 deficiency of the Gram-Charlier distribution in 
SUr<£L zone applications. 

2. Skewness and Kurtosis 

The physical geometry of the waves, which 

characteristically have narrow, steep crests and wide, 


Wemovwecreuchs ain the surf zone produce a positive 


ols 








skewness, indicating a greater amount of time below, 
rather than above, the mean water level. This is true 
except at Wave Gauge #3 on 6 March. (See Table II and 
Appendix F.) Data from this wave gauge was analyzed on 
this date only, because visual observations noted that 
most waves were actually breaking on this gauge, rather 
than on Wave Gauges #1 and #2 as was otherwise the case. 

It should follow that the wave induced particle 
velocities would have probability density functions similar 
to the waves. However, of the eleven horizontal velocity 
data records, four did not have positive skewness. These 
deviations occurred at Flow Meter #1 on 4 March and 
29 May and Flow Meter #2 on 5 and 6 March. Reflected 
Waves were originally thought to be a feasible explanation 
for the negative skewness. If this were so, Flow Meter #2 
on 4 March and 29 May would also expect to have negative 
Values. This was not the case. The reason for the large 
skewness and kurtosis values for Wave Gauge #1 on 5, 6 and 
8 March is not known. Similar anomalous results have been 
found in deeper water outside the surf zone by Thornton 
poe krapohl (1974). 

Kurtosis indicates the peakedness of each 
Meameter. Just prior to breaking the waves achieve the 
greatest degree of peakedness and should have the highest 
hur tosis value at this time. Visual observations noted 
mat mOSt waves were breaking at Wave Gauge #1 on 


@ March and 29 May, at Wave Gauge #2 on 5 March and on 


oy 





Wave Gauge #3 on 6 March. It should be expected that the 
kurtosis values would be greatest at these given gauges 
aoamecne respective dates. With the exception of 29 May, 
mis did not turn out to be true. The wide surf zone 
caused by varying breakers heights on 8 and 9 March make 


poeration Of kurtosis values difficult on these dates. 


fee COMPARISONS OF THEORETICAL AND MEASURED POWER SPECTRA 
The theoretical velocity spectra generated from the 
wave spectra and the measured velocity spectra were 
Sened if Order to determine the applicability of linear 
feieey) wave theory in the breaker zone. Thornton and 
Richardson (1974) measured waves with maximum heights of 
meomeend Showed that the application of a linear theory 
meatener £Unction to the surface profile spectrum resulted 
mecca l culation of the horizontal water particle velocity 
that underpredicted measured wave-induced horizontal 
femeerty spectral components by about 50 percent. The 
coherence values between waves and horizontal velocity 
mememnigh, ranging above 0.75 for an approximate frequency 
mamee Of 0.075 to 0.60 Hz. The phase angle computations 
Showed the calculated velocity components leading the 
measured velocity components by an average of 20 degrees. 
Bub (1974) took meaSurements in mild surf resulting in 
fmteoretical horizontal velocity values 13 percent lower 
@ean measured velocity. Coherence values varied from 
(moto nearly 1.0 over a frequency range of 0.10 to 0.65 Hz. 


An approximate zero degree phase was observed out to 1.0 Hz. 
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For this study wave heights of approximately 1 m 
were typically observed at Del Monte Beach and up to 2 m 
were recorded at Carmel River Beach. The theoretical 
horizontal velocity spectra were computed from the wave 


Spectra using the linear theory transfer function 


Pome esh kK (htz))- 2 
ery SD) (2) 
where 
ye = calculated horizontal velocity power 
spectrum (m*/sec), 
Sy (4) = measured wave height power spectrum (m*-sec), 


G6 = angular frequency = anf (sec +), 
k = wave number 2). 
h = mean water depth (m), 
z = depth of flow meter below the mean water 
depth (m). 
Meemnes 8 and 9 are graphs of the measured and calculated 
Spectral velocity values at Tower #1 for 4 March and 
29 May. Additional Spectra are contained in Appendix G. 
On the average at Del Monte Beach, the calculated hori- 
moneal velocity values underestimated the measured veloci- 
ties by 79%. The same theoretical calculations for the 
data collected at Carmel River Beach underestimated the 
measured values by 86%. 
A possible explanation for the low calculated values 


is that the wave profiles change during the shoaling and 
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igeakine process, DeECcCOMming more peaked at the crest, 
faseilessening the applicability of linear (Airy) theory 
which assumes sinusoidal wave shapes. Additionally, the 
wider separation between the measured and theoretical 
mpeetral values above 0.5 Hz may be explained by turbu- 


lence of the breaking waves in the higher frequency areas. 


eee XAMINATION OF SPECTRAL HARMONICS 

An examination of the low frequency spectral components 
for each data set indicates that pronounced approximate 
Narmonics are present in the data analyzed from Del Monte 
Beach. Such harmonics are not evident in the Carmel River 
Peace spectra. The frequencies and amplitudes of each 
fimeryzed Spectral peak 1s listed in Table III. The fre- 
Miiemey ditterences between Spectral peaks and the ratios 
gumenme Second, third and fourth peaks versus the primary 
Meameare listed in Table IV. 

Spectral analysis shows that definite second, third 
and fourth harmonics are present for 4 March and 5 March 
with few exceptions. Evidence of second and third har- 
Meomics 1s present in the analyzed data from Wage Gauges 
feamd #2 On 6, 8 and 9 March. An examination of the 
eieectral amplitudes shows that there is no consistent 
relationship between various components either for a 
given data set or when comparing like components of 
feeterent data sets. 

The maximum wave energies were concentrated at frequen- 


Siccmraneime iron 0.061 Hz to 0.092 Hz corresponding to 


a) 
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Peetods Of 16.4 to 10.9 seconds, respectively, which agreed 
with visual observations. Additionally, prominent sub- 
harmonic peaks between 0.011 and 0.022 Hz are noted in the 
measured velocity spectra computed from Del Monte Beach 
tere this agrees with the work of Bub (1974), who 
observed a sub-harmonic at 0.011 Hz. Guzza and Davis (1974) 
State that in a theoretical analysis of edge waves excited 
by incoming waves, the prominent edge wave mode is the 

first sub-harmonic of the primary frequency. It is plausi- 
Ble, therefore, that the observed low frequency sub-harmonic 
peaks were caused by edge waves in the surf zone. Such 
sub-harmonics are not observed in the measured velocity 


spectra generated from Carmel River Beach data. 


Eeeeeaase AND COHERENCE 

The curling crests of the unstable breaking waves 
Bemgmuerlcad the horizontal particle velocities in the body 
of the waves as can be seen in the phase measurements 
between waves and velocities. (See Table V.) The maximum 
phase difference would be expected at breaking. Data from 
@ower #1 0n 4 March indicates that the theoretical hori- 
zontal velocity leads the measured horizontal velocity by 
about 20 degrees. Since most waves observed were breaking 
on Tower #1, a possible explanation for the 30 degree phase 
angle at Tower 2 is that the waves reorganized between the 
towers and broke a second time at Tower #2. Since this 


author's observations concentrated primarily on the initial 
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breaking position, the possible reorganization and second 
breaking were not observed. The phase angles on 5 March 
Gee and 20 degrees at Towers #1 and #2, respectively. 
This agrees with observations and indicates that the 
breakers started to curl at Tower #1 and broke at Tower #2. 
Wiemzero degree phase relationships on 6 March show that 
the waves had not started to break until shoreward of 
Tower #2. This agrees with the observation of waves 
breaking on Tower #3. The phase angles on 8 and 9 March 
indicate that the smaller breakers began to curl at Tower #1 
and were near breaking at Tower #2. Again, this agrees 
Wap ODServations. 

iacontrast to the plunging-type breakers at Del Monte 
Beach, surging-spilling breakers were observed at Carmel 
Peegeme beach, This was to be expected considering the beach 
Eeomes. A spilling breaker is characterized by turbulent 
water forming at the wave crest and eventually flowing down 
and covering the leading edge of the wave. Spilling begins 
at the crest when a small tongue of water moves faster 
than the wave form as a whole (Galvin, 1972). The measured 
phase angle at Carmel River Beach shows that the cascading 
@rests led the horizontal particle velocities in the body 
of the wave by 10 degrees. 

The values of coherence show varying degrees of linearity 
mmetne transter functions between different spectra. Co- 
ference between wave and horizontal velocity spectra at 


Tower #1 on 29 May is above 0.53 out to 0.58 Hz (See 
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lpeore 9), above 0.68 out to 0.53 Hz on 8 March and above 
wou out to 0.41 Hz on 9 March. (See Appendix G.) Other 
coherence calculations show a fair to poor degree of 
Mmenearity, which is consistent with the observed strong 


non-linearities. 


feeeecORETICAL AND CALCULATED PHASE SPECTRA 

A theoretical phase spectrum was calculated for com- 
parison with the measured wave spectrum from two wave 
gauges. The theoretical phase was computed from the 
oma ttonship 


oo (3) 


where 
od = phase, 
Oo = angular frequency = 2rf Gaes- 
x = distance between wave gauges (m), 


1/2 (m/sec), 


Peemecclenity — (chi) 

g = gravity hilsee 

h = mean water depth (mn). 
Figure 10 shows typical examples of the theoretical and 
measured phases. The theoretical phases computed from Del 
Monte Beach data underpredict the actual phase relationships, 
whereas an overprediction is noted in the Carmel River Beach 
phases. The author then determined what value of celerity 


(C,) would cause the theoretical and actual phases to 


agree out to the first zero crossing (where $¢ = 180° = 7) 
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Dyaetirst solving for C, poten iteration Of equation (3) 


C, = 21x (4) 


and then solving 


C, = (g(hta))1/? (5) 


meme. lhe results are shown in Table VI. No evident 
Pelartronshiip between a and other parameters, such as beach 
Slope, mean water depth and standard deviation, could be 


determined. 


45 





SNOSTUVdNOO dSVHd YOd CALVTNOTVO SYALANVAVd 





Selle: Gaeetoliciey ls 





iy ae ees LOG Va Sis T8829 °T ey C# OM T# OM ABW 6Z 
0998°0 pso'V SO9LT*O Co et pssvs’o0 bas. Os C# OM - T# OM YteWN 6 
Deo 0 688 °¢ Sees be (0 09°TT 69906°0 to 2c C# OM T# OM YoteN 8 
S56 0720 Ss ae SvOOT’O ONG ane CO6O0T’T tee ES C# OM - 2# OM YteW 9 
9880°0 p8s°¢ Eyes 0 OC iT 890SP°T T:8°SZ C#H OM T# OM YrIteW 9 
Get 0 pSs'¢ veo 0 09°TT SZ60b'T ieee el C# OM T# OM UYrYteW S 
S060 Soe) tS SZ79ST°O CT HE O20. 0) tee vi C# OM T# OM YSILeNW PV 
(Si9z9uW) 0 (D905 /u) 9 Cables (Si9}o9W)x (Ssiozouw) yY OS SJUSUNIISUT 518q 


Yovog 


406 





V. CONCLUSIONS 


Breaking waves in the surf zone can be characterized 
as being highly non-linear. Qualitative observations of 
wave and velocity profiles show secondary waves indicative 
Of non-linear waves. The probability density functions 
calculated for the flow velocity records compared better 
Peeeaeetice Gram-Charlier distribution than with the Gaussian 
distribution when tested using the chi-square goodness-of- 
fit test. However, seven of the thirteen wave records 
more closely approximated the Gaussian distribution. This 
[Pmienesult of the truncated form of the Gram-Charlier 
Pe bability density function used in calculations, and 
Memes Out the importance of including higher order moments 
in describing wave phenomena in very shallow water. 

ities vaiues of the horizontal power spectral components 
calculated from wave spectra using linear theory indicate 
peegetitative, but not a quantitative, relationship. Com- 
paring the results of Bub (1974), Thornton and Richardson 
(1974) and this author, it 1s evident that increased 
breaker heights result in larger estimation errors. 
Additional investigation is warranted to determine if 
specific breaker heights and/or types result in specific 
M@erestimation values. If so, linear theory may be 
applied to achieve a reasonably accurate determination of 


flow velocities. The coherence values between waves and 
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bememzontcal vellocity decrease as breaker height increases, 
indicating that the wave motion becomes more non-linear 
and more turbulent under this condition. 

The theoretical phase calculations based on the linear 
theory approximation for celerity did not accurately 
predict the observed phases in this research. The values 
of a, determined so as to force the theoretical phases to 
agree with the observed values at the first zero crossing, 
eimow moO consistency. 

Use of electromagnetic flow meters and the improved 
Capacitance wave gauges and instrument towers permits the 
gathering of accurate continuous data with sturdy and 


reliable equipment. 
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APPENDIX B 
CAPACITANCE WAVE GAUGE CALIBRATIONS 
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APPENDIX D 


ANALYSIS DETAILS AND FLOW CHART 


A complete flow chart of data acquisition and analysis 
is presented on the following page. The parameters used 
meemme spectral analysis were chosen to balance the longest 
record that the computer could reasonably analyze with the 
best resolution over the frequency range of interest. 
Additionally, consideration was given to the computer run 
ELMe. 

Im order to determine the optimum Nyquist frequency for 


data analysis the following equations were used: 
At = (VDT) (NSKIP) (NCHAN) (ley 


and 


f= XU (2) 
where 
At = sampling interval, 
VDT = sampling rate, 
NSKIP = number of samples skipped in 
ot Gade dawdeeat fay , 
NCHAN = number of channels, 
fy = Nyquist frequency. 
By varying the parameters of equation (1) various Nyquist 
frequencies were obtained and used during preliminary IBM 


500 computer analysis. 
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CALIBRATION FACTORS 


Date Instrument ae 


(meters) (meters/volt) 

4 March WG #1 =o ees 
and WG #2 =128.70 4.19 

5S March WG #3 -1.800 Paes 
© March WG #1 = bea 2 Paes 
WG #2 2412020 4.19 

WG #3 -1 .800 y (ae be 

8 March WG #1 = eso ees 
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WG #3 -1.810 i Boe 

Sievancch WG #1 -1.492 LZ 
WG #2 Ss 0: 4.19 

WG #3 -1.830 1.14 

29 May WG #1 -0.028 0.69 
WG #2 ~0.370 0.66 

All dates FM #1 10 52 uUS 
FM #2 0.0 Sail 
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POWER, COHERENCE AND PHASE SPECTRA 
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